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Abstract We review some recent studies on the string model of confinement inspired 
by the strong-coupling regime of QCD and its application to exotic multiquark config- 
urations. This includes two quarks and two antiquarks, four quarks and one antiquark, 
six quarks, and three quarks and three antiquarks with a careful treatment of the 
corresponding few-body problem. 
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The physics of multiquark systems remains rather delicate. On the experimental 
side, most of the candidates which have been announced are not yet confirmed [T]. 
Moreover, the search has been focused mainly on states with hidden flavor, which might 
be either genuine multiquarks or ordinary hadrons strongly coupled to their decay 
channels. Other sectors such as double charm mesons have not been much explored. 

As for the theory, the main question is how to extrapolate quark models adapted to 
mesons and baryons toward the multiquark sector. Then it remains to solve carefully the 
few-quark problem. Many configurations lie at the edge between binding and instability, 
so the wave function should incorporate both short-range components where the quark 
interaction takes place and long-range components which include the clusterized part 
of the wave function. 

In the studies summarized here, we focused on the role of confinement, with a 
string model inspired by [2,3 . The linear potential corresponds to the minimal length 
for the gluon flux tube linking the quark to the antiquarks. For baryons, it becomes 
the y-shape interaction which minimizes the cumulated length of the flux tubes. For 
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Table 1 Four-quark energies, E4, for the different confinement models described in Eq. |3j 
compared to their threshold, T4, as a function of the mass ratio M/m. 



M/m 




E 4 (QQqq) 




T 4 (QQqq) 


E' A {QQqq) 


Tl(QQqq) 




Vf 


v b 


V s 




Vf 




1 


4.644 


5.886 


4.639 


4.676 


4.644 


4.676 


2 


4.211 


5.300 


4.206 


4.248 


4.313 


4.194 


3 


4.037 


5.031 


4.032 


4.086 


4.193 


3.959 


4 


3.941 


4.868 


3.936 


3.998 


4.117 


3.811 


5 


3.880 


4.754 


3.873 


3.942 


4.060 


3.705 



tetraquarks and beyond, the potential is mainly of the flip-flop terms with an inter- 
action inside each cluster and a minimization over all possible cluster decompositions. 
There are connected Steiner diagrams which link all quarks and antiquarks, see Fig. [I] 

In Refs. [4,5,6 , we solved the tetraquark, pentaquark and hexaquark problems with 
this confining potential alone, for systems whose quarks are assumed to be different, 
even in the case where they bear the same mass. This means that are neglected alto- 
gether antisymmetrization, relativistic effects, short-range corrections, spin-dependent 
effects, etc. We are now studying the interplay with chromomagnetic forces and the 
role of antisymmetrization. 

This string model replaces the "color-additive" ansatz that was widely used in the 
early days of hadron spectroscopy and reads 

V(r 1 ,...) = -^^2~X i :\ j v(r ij ) . (1) 

i<j 

For baryons, the color-additive model gives the well-known "1/2" rule between the 
quark-quark and the quark-ant iquark potentials [7 , while the string model leads to [8 

y(n,r 2 ,r 3 ) = Amin^! + r k2 + r k3 ) , (2) 

k 

reminiscent of the celebrated problem of Fermat and Torricelli, later generalized as the 
Steiner problem [9]. 

For tetraquarks, the potential is the minimum of the flip-flop and connected Steiner- 
tree terms. The former is the minimum of the two possible configurations made of two 
mesons, and the later is a generalization of the y-shape interaction, with optimization 
with respect to two junction points. More precisely, 

V s =mm(V f ,V b ) , 

(3) 

Vf = Amin(ri 3 + r 2 4,r 2 3 + n 4 ), V b = Amin(r lfc + r 2 k + r ki + r i3 + r i4 ) . v ; 

k,£ 

For a detailed discussion see Ref. [4]. 

The results are shown in Table [l] Clearly, as M/m increases, a deeper binding is 
obtained for the flavor-exotic (QQqq) system. For the hidden-flavor (QQqq), however, 
the stability deteriorates, and with our variational approximation, for M/m > 1.2, the 
system becomes unbound with respect to its lowest threshold (QQ) + (qq)- 

The pentaquark potential of Ref. [5] is similarly a combination of flip-flop and 
connected terms 



V P = min(y ff ,y st ) , V s = mm[r li + V Y (r j ,r k ,r i )] , 



(4) 
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Table 2 (Q 3 q 3 ) and (Q 3 q 3 ) variational energies E, compared to their threshold energy T, 
as a function of the mass ratio M/m. A = E — T is the energy difference. 



M/m 


E(Q 3 q 3 ) 


T(QV) 


A(Q 3 q 3 ) 


E(Q 3 q 3 ) 


T(Q 3 q 3 ) 


A(Q 3 q 3 ) 


1 


7.237 


7.728 


-0.491 


6.981 


6.981 


+0.000 


2 


6.524 


6.929 


-0.405 


6.314 


6.335 


-0.021 


3 


6.209 


6.543 


-0.334 


6.030 


6.079 


-0.049 


4 


6.014 


6.298 


-0.294 


5.868 


5.940 


-0.072 


5 


5.890 


6.123 


-0.233 


5.762 


5.852 


-0.090 



with {z,j, k,£} suitably permuted, and V$t denoting the connected Steiner trees. 

The calculation was performed using a hyperspherical approximation and a Gaus- 
sian expansion. Both methods agree that systems made of four light quarks and one 
light antiquark are clearly below the lowest meson-baryon threshold. This is the oppo- 
site to what occurs when a pair- wise interaction is considered. 

If the calculation is repeated in the case of four unit masses and an infinitely massive 
antiquark or quark, the pairwise model gives an energy well above the threshold for both 
(Qqqqq) and (qqqqQ). On the other hand, the flip-flop potential gives E^(Qqqqq) — 
E t h(Qqqqq) ~ —0.114, and E^(qqqqQ) — E^qqqqQ) ~ —0.2777, which is sufficient 
to demonstrate their stability. 

For the (q 6 ) configurations, there are again two types of diagrams: flip-flop and 
connected Steiner trees. The potential is minimized with respect to all permutations. 
For the (q 3 q 3 ) states, there are several possibilities: flip-flop with either a baryon and 
an antibaryon, or three mesons, or a meson and a tetraquark, and also some connected 
diagrams with four or more junctions. The dynamics is dominated by the flip-flop terms, 
while the connected diagrams with y-shape junctions play a minor role for binding. 
Moreover, the dynamics of baryon is qualitatively similar with a pair- wise potential 
rjj/2 and the y-shape model, this leading to easier computations without spoiling 
the spirit of the model. For a more detailed discussion see Ref. [BJ. 

The results are shown in Table[2] For the "dibaryon", (Q 3 q 3 ), the system is found 
stable against dissociation into two baryons. However, the stability deteriorates when 
the mass ratio increases. The behavior is reasonably linear and therefore the limit where 
the system becomes unbound can be estimated to be of the order of M/m ^8 — 10. 
Hence a triple-charm dibaryon is predicted but not a triple-beauty one. 

In the case of (Q 3 q 3 ), this is more delicate. For M = m = 1, the results in 
Table [2] suggest there is a shallow bound state below the lowest threshold. This means 
that the effective interaction between the (qq) mesons is attractive. Now, as M/m 
further increases, the (Q 3 ) + (q 3 ) threshold will become degenerate with the lowest 
mesonic threshold. This will favor binding, as the six-body wave function will contain 
two different decompositions into clusters with relative motion that will interfere to 
improve binding. However, for even larger values of the mass ratio M/m, no multiquark 
configuration can acquire enough binding to compete with the compact (Q 3 ), and 
the system becomes unstable against rearrangement into (Q 3 ) + (q 3 )- Perhaps, some 
metast ability could be observed with respect to some higher threshold. 

Let us summarize and suggest some possible further studies. 

1. The string model of confinement which combines flip-flop and connected flux 
tubes of minimal length gives more attraction than the additive pairwise model with 
color factors that was used in early multiquark calculations. 

2. This potential is flavor independent. By changing the constituent masses in the 
kinetic-energy part of the hamiltonian, one can modify the binding. For tetraquarks, 
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Fig. 1 A few examples of flip-flop terms and connected Steiner trees. From left to right, top 
to bottom: the three contributions to the tetraquark, the flip-flop term of the pentaquark (to 
be minimized by permuting the quarks), and some disconnected and connected contributions 
to the (Q 3 q 3 ) confinement. 

(QQqq) is more stable when the mass ratio M/m increases, while (QQqq) becomes 
unstable. For pentaquarks the calculation has been restricted to five equal and finite 
masses, or one infinitely massive constituent surrounded by four equal masses. The 
property of stability was found to survive for both (Qqqqq) and (qqqqQ). In the six- 
quark case the binding energy decreases for (Q 3 g 3 ), while for (Q 3 g 3 ), it first increases 
and then decreases. In the limit of large M/m, no six-body configuration can compete 
with the deep binding of (Q 3 ) that enters the lowest threshold. 

3. These potentials can be seen as a Born-Oppenheimer limit. Our aim is to refor- 
mulate the interaction as an operator in color space, to go beyond the adiabatic limit 
and include the constraints of antisymmetrization. 
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